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Abstract 

Background: Mycobocterium bovis is the causative agent of bovine tuberculosis (BTB), a pathological infection with 
significant economic impact. Recent studies have highlighted the role of functional genomics to better understand 
the molecular mechanisms governing the host immune response to M. bovis infection. Furthermore, these studies 
may enable the identification of novel transcriptional markers of BTB that can augment current diagnostic tests 
and surveillance programmes. In the present study, we have analysed the transcriptome of peripheral blood 
leukocytes (PBL) from eight M. bows-infected and eight control non-infected age-matched and sex-matched 
Holstein-Friesian cattle using the Affymetrix® GeneChip® Bovine Genome Array with 24,072 gene probe sets 
representing more than 23,000 gene transcripts. 

Results: Control and infected animals had similar mean white blood cell counts. However, the mean number of 
lymphocytes was significantly increased in the infected group relative to the control group {P = 0.001), while the 
mean number of monocytes was significantly decreased in the BTB group {P = 0.002). Hierarchical clustering 
analysis using gene expression data from all 5,388 detectable mRNA transcripts unambiguously partitioned the 
animals according to their disease status. In total, 2,960 gene transcripts were differentially expressed (DE) between 
the infected and control animal groups (adjusted P-value threshold < 0.05); with the number of gene transcripts 
showing decreased relative expression (1,563) exceeding those displaying increased relative expression (1,397). 
Systems analysis using the Ingenuity® Systems Pathway Analysis (IPA) Knowledge Base revealed an over- 
representation of DE genes involved in the immune response functional category. More specifically, 64.5% of genes 
in the affects immune response subcategory displayed decreased relative expression levels in the infected animals 
compared to the control group. 

Conclusions: This study demonstrates that genome-wide transcriptional profiling of PBL can distinguish active M 
5o\//5-infected animals from control non-infected animals. Furthermore, the results obtained support previous 
investigations demonstrating that mycobacterial infection is associated with host transcriptional suppression. These 
data support the use of transcriptomic technologies to enable the identification of robust, reliable transcriptional 
markers of active M. bovis infection. 
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Background 

Bovine tuberculosis (BTB) poses a serious threat to the 
health of domestic cattle herds worldwide. Infection is 
caused by the bacterium Mycobacterium bovis, an intra- 
cellular pathogen closely related to Mycobacterium 
tuberculosiS'the causative agent of human tuberculosis. 
M. bovis infection is often slow and progressive with 
Umited cUnical symptoms. Although improved diagnos- 
tic tests and slaughter policies have done much to con- 
trol and reduce the incidence of infection, BTB has 
remained recalcitrant to eradication in many countries 
where control programmes have been implemented 
[1-3]. 

Failure to detect and remove all infected animals from 
herds is partly due to limitations in the sensitivity of the 
current diagnostic tests, which often comprise an in vivo 
single intradermal comparative tuberculin test (SICTT) 
performed alone, or in combination with an in vitro 
enzyme-linked immunosorbent assay (ELISA)-based test 
for interferon gamma (IFN-Y)-an established biomarker 
of mycobacterial infection [4-6]. Diagnoses can be 
further confounded by exposure to environmental non- 
pathogenic mycobacterial antigens, which can generate 
false SICTT-positive signals in cattle [7]. Protection 
from natural M. bovis infection in cattle may be 
achieved through vaccination with M. bovis bacillus 
Calmette-Guerin (BCG); however, the level of protection 
attained is variable. In addition, current diagnostics can- 
not effectively differentiate between M. bovis -infected 
and BCG-vaccinated animals, thus compromising man- 
agement strategies [8]. Consequently, there is a pressing 
need for novel M. bovis diagnostic methods with 
increased sensitivity and specificity. 

The host immune response to mycobacterial infection 
is a complex process that involves interaction between 
the innate and adaptive immune systems. Upon initial 
exposure (generally via inhalation), bacilli are phagocy- 
tosed by host alveolar macrophages, which recognise 
mycobacteria using a diverse range of pathogen recogni- 
tion receptors (PRRs), such as the Toll-like receptors 
(TLRs) and the nucleotide-binding oligomerisation 
domain (NOD)-like receptors (NLRs) [9-13]. Activation 
of macrophage PRR-mediated signalling pathways result 
in the release of endogenous cytokines, which initiate an 
adaptive immune response characterised by the secre- 
tion of proinflammatory cytokines, such as IFN-y and 
tumour necrosis factor (TNF-a), by activated T cells 
[14]. In particular, IFN-y activates infected macrophages 
and enables the formation of granulomas-coUections of 
inflammatory cells comprising T cells, B cells, non- 
infected macrophages and neutrophils, which surround 
infected macrophages and act as barriers to contain and 
prevent dissemination of the infection [15]. In most 
cases, the host innate and adaptive immune systems 



successfully control mycobacterial growth within granu- 
lomas resulting in asymptomatic latent infection [16,17]. 
However, in some cases impairment of immune function 
can result in the development of active tuberculosis 
leading to disease progression [3,17-20]. 

Recently, functional genomic technologies have been 
used to investigate the molecular mechanisms and cellu- 
lar pathways underlying the host immune response to 
mycobacterial infection [[21,22], for reviews see [23,24]]. 
Furthermore, results from these studies have the poten- 
tial to identify molecules that are critical for host/patho- 
gen survival during infection, and which may serve as 
robust, reliable transcriptional markers of mycobacterial 
infection [22]. 

Previously, we investigated the transcriptional profiles 
of peripheral blood mononuclear cells (PBMC) from M. 
Z70v/5-infected and non-infected control animals using 
the immuno-specific BOTL-5 microarray (containing 
1,391 gene probe sets; Gene Expression Omnibus [GEO] 
accession: GPL5751) and showed that suppression of 
innate immune genes was associated with BTB [25]. In 
the current study, we extend this earlier work by investi- 
gating the transcriptional profile of peripheral blood leu- 
kocytes (PBL) isolated from eight M. ^ov/5-infected and 
eight non-infected control animals using the genome- 
wide high-density Affymetrix® GeneChip® Bovine Gen- 
ome Array. These 16 animals were sampled specifically 
for the present study and have not been used for any 
previous research work. The Affymetrix® GeneChip® 
Bovine Genome Array contains 24,072 gene probe sets 
representing more than 23,000 gene transcripts http:// 
www.affymetrix.com. In addition, we have adopted a 
systems biology approach using the Ingenuity® Systems 
Pathway Analysis (IPA) Knowledge Base http://www. 
ingenuity.com for analysis of both over-represented cel- 
lular functions and known molecular canonical pathways 
from the resulting gene expression data. 

The results presented in the current study contribute 
a novel layer of information regarding the gene expres- 
sion profile of PBL from M. bovis-infected animals and 
highlight the value of high-throughput genomic technol- 
ogies in understanding the host immune response to 
BTB. Furthermore, these results may facilitate the devel- 
opment of novel diagnostics for the detection of M. 
bovis infection in domestic herds. 

Methods 

Experimental animals 

Sixteen age-matched female Holstein-Friesian animals 
from cattle herds that had not been analysed previously 
were used for this study. Eight infected individuals were 
selected from a panel of naturally M. bovis -infected ani- 
mals maintained for on-going disease surveillance at the 
Irish Department of Agriculture, Fisheries and Food, 
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Backweston Laboratory Campus (Celbridge, Co. Kildare, 
Ireland). These animals had a positive single intradermal 
comparative tuberculin test (SICTT) result where the 
skin-fold thickness response to purified protein deriva- 
tive (PPD)-bovine exceeded that of PPD-avian by at 
least 12 mm. All of these animals were also positive for 
the whole blood IFN-y-based BoviGAM® assay [Prionics 
AG, Switzerland] (data not shown). In addition, these 
cattle were confirmed for BTB following detailed post- 
mortem pathological examination and/or culture. 
Briefly, bronchial, mediastinal, submandibular, retro- 
pharyngeal, mesenteric and hepatic lymph nodes and 
lungs were examined macroscopically for tuberculosis 
lesions. Suspected lesions were cultured on Stonebrinks 
and Lowenstein- Jensen media at 37°C for eight weeks to 
detect M. bovis [26] . Non-infected control animals were 
selected from a herd with no recent history of M. bovis 
infection. The control animals were shown to be nega- 
tive for both the SICTT and IFN-y tests (data not 
shown). All animal procedures detailed in this study 
were carried out according to the provisions of the 
Cruelty to Animals Act (licenses issued by the Depart- 
ment of Health and Children) and ethics approval for 
the study was obtained from the UCD Animal Ethics 
Committee. 

Blood collection 

Two 8 ml vacutainers® (Becton-Dickinson Ltd., Dublin, 
Ireland) of heparinised blood were collected from each 
animal, approximately 12 months after positive SICTT 
testing. One vacutainer® was retained for haematological 
analysis using a Cell-Dyn 3500 haematology analyser 
(Abbott Laboratories Ireland Ltd., Dublin, Ireland); all 
haematological analysis was performed using 1 ml of 
blood. The other vacutainer® was used for RNA isola- 
tion from peripheral blood leukocytes (PBL); the whole 
white blood cell fraction consisting of T and B lympho- 
cytes, NK cells, monocytes, neutrophils, basophils and 
eosinophils. The count data from the leukocyte cell 
populations of infected and non-infected animals were 
assessed using the two-sample, two-tailed Student's t- 
test, following Kolmogorov-Smirnov tests of normality 
and Levene's F-test for equality of variance using the 
Minitab statistical package version 16 (Minitab Ltd., 
Coventry, UK). 

RNA extraction and microarray analysis 

All RNA extractions were performed within two hours 
of blood collection. Briefly, 7.5 ml of whole heparinised 
blood was mixed with 42.5 ml of erythrocyte blood lysis 
buffer (10 mM KHCO3, 150 mM NH4CI, 1 mM EDTA 
pH 8.0), and incubated for 5 min at room temperature 
with gentle agitation. Following centrifugation (750 g for 
10 min) the pelleted cells were washed once with Ix 



phosphate buffered saline (Invitrogen Ltd., Paisley, UK). 
The cell pellet was then fully resuspended in 2 ml Tri- 
zol® reagent (Invitrogen Ltd., Paisley, UK) and RNA was 
extracted as per the manufacturer's instructions. The 
RNA was further purified using an RNeasy® kit with 
on-column DNase treatment (Qiagen Ltd., Crawley, UK) 
according to the manufacturer's instructions. RNA 
quantity and quality was assessed using both the Nano- 
j)j.QpTM 1000 spectrophotometer (Thermo Fisher Scien- 
tific, Waltham, MA, USA) and the Agilent 2100 
Bioanalyzer using an RNA 6000 Nano LabChip kit (Agi- 
lent Technologies, Cork, Ireland). All samples displayed 
a 260/280 ratio greater than 1.8 and RNA integrity 
numbers (RIN) greater than 8.0. 

cDNA labelling, hybridisation and scanning for the 
microarray experiments were performed by Almac Diag- 
nostics (Craigavon, Co. Armagh, Northern Ireland) 
using a one-cycle amplification/labelling protocol on the 
Affymetrix® GeneChip® Bovine Genome Array (Affyme- 
trix UK Ltd., High Wycombe, UK). 

Statistical analysis of microarray data 

Affymetrix® GeneChip® Bovine Genome Array data 
were analysed using Bioconductor [[27]; http://www. 
bioconductor.org] contained within the R statistical 
package http://www.r-project.org. Normalisation of raw 
data was performed using the Factor Analysis for 
Robust Microarray Summarization (FARMS) algorithm. 
The FARMS algorithm uses only perfect match (PM) 
probes and a quantile normaUzation procedure, provid- 
ing both P- values and signal intensities [28]. In addi- 
tion, gene expression profiles for each animal were 
clustered using the Hierarchical Ordered Partitioning 
and Collapsing Hybrid (HOPACH) clustering algo- 
rithm in Bioconductor with Euclidean distance as the 
distance metric [29]. 

Normalised data were then further subjected to filter- 
ing for informative probes sets using the R package II 
Nl-calls [30]. This defines a probe set as being informa- 
tive when many of its probes reflect the same change in 
mRNA concentration across arrays. Differentially 
expressed genes were extracted using the Linear Models 
for Microarray Data (LIMMA) package http://biocon- 
ductor.org/packages/release/bioc/html/limma.html con- 
tained within the R statistical package. Genes displaying 
differential expression patterns between control and 
infected groups were annotated using the Affymetrix® 
bovine gene annotation http://www.affymetrix.com. The 
Benjamini-Hochberg multiple -testing correction method 
[31] was applied to all differentially expressed genes to 
minimise the false discovery rate (FDR) and adjusted P- 
values for differentially expressed (DE) genes were cal- 
culated. For genes represented by multiple probe sets 
the mean expression value is reported. Bootstrapping for 
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the hierarchical clustering analysis was performed using 
a custom Perl script (see Additional file 1). 

Systems biology analyses 

Ingenuity® Systems Pathway Analysis (IPA, Ingenuity 
Systems, Redwood City, CA, USA; http://www.ingenuity. 
com) was used to identify canonical pathways and func- 
tional processes of biological importance within the list 
of DE genes. The Ingenuity® Knowledge Base contains 
the largest database of manually-curated and experimen- 
tally-validated physical, transcriptional and enzymatic 
molecular interactions. Furthermore, each interaction in 
the Ingenuity® Knowledge Base is supported by pre- 
viously published information. 

For the IPA analyses, the Affymetrix® GeneChip® 
Bovine Genome Array was used as a reference gene set. 
All DE genes with an adjusted P value < 0.05 were 
included. For duplicate probe IDs, the average log2 
expression fold change was used. Only DE genes map- 
ping to molecules in the Ingenuity® Knowledge Base 
were used for systems analysis. Functional analysis of 
genes was performed using IPA to characterise the bio- 
logical functions of the DE genes between the BTB and 
control groups. For this, IPA performed an over-repre- 
sentation analysis that categorises the DE genes within 
the uploaded list into functional groups using the Inge- 
nuity® Knowledge Base. Each category in IPA is ranked 
based on the number of DE genes falling into each func- 
tional group. Right-tailed Fisher s exact tests were used 
to calculate a P-value for each of the biological function 
assigned to list of DE genes. 

IPA contains a large library of known canonical path- 
ways that were overlaid with the DE genes to identify 
major biological pathways associated with M. bovis 
infection in PBL. The significance of the association 
between DE genes and the canonical pathway was 
assessed using two methods: (1) a ratio of the number 
of molecules from the DE gene data set that map to the 
pathway, compared to the total number of molecules 
that map to the canonical pathway based on the refer- 
ence gene list; and (2) a Fisher's exact test that generates 
a P-value for the assignment of the DE genes to a parti- 
cular canonical pathway compared to the reference gene 
list. Canonical pathways were then overlaid with the 
expression values of the DE genes. 

Real time quantitative reverse transcription PGR (qRT- 
PCR) validation of microarray results 

cDNA was prepared using 500 ng of total RNA from 
each sample from the microarray study using the High 
Capacity cDNA Reverse Transcription Kit (Applied Bio- 
systems, Warrington, UK) in a 20 (il reaction using ran- 
dom primers according to the manufacturer's 
instructions. cDNA was diluted 1:50 and stored at -20°C 



prior to performing real time quantitative reverse tran- 
scription PGR (qRT-PGR). 

Real time qRT-PGR reactions were performed using 
Fast SYBR® Green Master mix (Applied Biosystems, 
Warrington, UK) on a 7500 Fast Real-Time PGR System 
apparatus (Applied Biosystems, Warrington, UK). Each 
reaction (20 contained 5 \i\ of the diluted cDNA 
(equivalent to 2.5 ng of total RNA) and 300 nM final 
concentration each of forward and reverse primer. Addi- 
tional file 2, Table SI provides a complete list of primer 
sequences and the target accession numbers for each 
real time qRT-PGR amplicon analysed. Real time qRT- 
PGR primers were designed using the PrimerS software 
[32,33] and where possible intron-spanning primers 
were selected (see Additional file 2, Table SI). Negative 
real time qRT-PGR controls and a six-point, four-fold 
dilution series from pooled cDNA from all animals were 
included on every real time qRT-PGR plate and indivi- 
dual PGR efficiencies were determined from the stan- 
dard curves using the qbase^^^^ software package [[34]; 
Biogazelle NV, Zwijnaarde, Belgium]. 

The PGR thermal cycling program consisted of one 
cycle at 50°G for 2 min, one cycle at 95°G for 10 min, 
followed by 40 cycles at 95°G for 15 s and 60°G for 1 
min. A dissociation step was included to confirm ampli- 
fication specificity and real time qRT-PGR products 
were analysed on a 2% agarose gel to confirm the pre- 
sence of a single discrete amplicon of the correct size. 
All real time qRT-PGR data was analysed using the qba- 
se^^^^ software package with efficiency correction and 
normalization was performed using two reference genes: 
the 60S ribosomal protein L19 gene {RPL19) and pepti- 
dylprolyl isomerase A (cyclophilin A) gene (PPIA). The 
two reference genes were selected using the geNorm 
algorithm in the qbase^^^^ package from a panel of 
eight genes tested (geNorm M > 0.15). 

The qbase^^^^ package generated a calibrated normal- 
ized relative quantity (GNRQ) of gene expression for 
each of the analysed samples. Log2 GNRQ values for 
both the control and the M. ^ov/5-infected animals were 
used for statistical analysis for all genes. One sample 
Kolmogorov-Smirnov tests, performed using the SPSS® 
version 18 software package (SPSS® Inc., Ghicago, IL, 
USA) were applied to the residuals of the log2 GNRQ 
values for each sample prior to statistical analysis to 
ensure the data conformed to a normal distribution— no 
significant departures from normality were observed for 
any of the genes analysed {P > 0.05). Two-tailed, two- 
sample Student ^-tests were used to assess differences 
between infected and control groups based on log2 
GNQR values using SPSS and the statistical package 
contained within Microsoft® Excel 2010 (Microsoft 
Gorp., Redmond, WA, USA). In addition, Levene's F- 
test was applied to the infected and control group log2 
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CNRQ values to assess equality of variance within the 
two groups using SPSS prior to performing the two- 
sample ^-tests; equality of variance was observed for 
each sample group. 

Geometric mean fold-changes in gene expression for 
the BTB group were calculated by dividing the geo- 
metric mean CNRQ value for the M. ^ov/5-infected 
group by the geometric mean CNRQ value for the con- 
trol group. The negative reciprocal in fold-change is 
reported where decreased gene expression was observed 
in the BTB group relative to the control group. 

Results 

Analysis of leukocyte populations in control and M. bows- 
infected blood samples 

The infected animals used in this study were chosen on 
the basis of their large responses to the comparative 
tuberculin skin test. The IFN-y levels measured in whole 
blood of the infected animals were at least 25-fold 
greater than in the healthy control cattle {P < 0.001, 
data not shown), demonstrating that the infected ani- 
mals were generating strong cell-mediated immune 
responses. At post-mortem, each of the infected animals 
displayed gross tuberculosis lesions in the lungs and 
thoracic lymph nodes and were classified as being in the 
advanced stage of clinical disease. 

To assess potential changes in leukocyte composition 
between control and M. bovis- infected samples, whole 



blood samples were subjected to haematological analy- 
sis (Figure 1). No significant difference was observed 
in the total white blood cell (WBC) count between 
control and M. bovis-'miQcted animals {P = 0.18). How- 
ever, significant increases in the mean number of lym- 
phocytes {P = 0.001) and significant decreases in the 
mean number of monocytes {P = 0.002) were observed 
in M. boviS'infected animals relative to the control 
animals. No significant differences in the mean number 
of eosinophils and neutrophils were observed between 
the two sample groups {P = 0.51 and P = 0.37, 
respectively). 

Summary of differentially expressed genes between 
control and M. bows-infected animals identified from 
Affymetrix® GeneChip® analysis 

The expression data generated for the current study are 
MIAME-compliant [35] and were deposited in the NCBI 
Gene Expression Omnibus (GEO) repository [36] with 
experiment series accession GSE33359. All array tran- 
scripts that passed informative probe filtering were used 
for cluster analysis to examine the grouping of samples 
based on infection status (Figure 2). In total, 5,388 tran- 
scripts passed the filtering process and were used for 
the cluster analysis, which showed a clear partitioning of 
samples based on their disease status, indicating a dis- 
tinct difference in expression profile between the two 
sample groups. The division of the two sample groups 
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Figure 1 Mean leukocyte cell population subset counts in control and M. 6ows-infected animals (n = 8 each group). Error bars represent 
the standard error for each mean cell count estimate. 
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Figure 2 Hierarchical clustering dendrogram generated using the complete gene expression profiles for the M. 6oWs-infected and 
control animals. The dendrogram was generated from the complete gene expression dataset for each animal using the Euclidean distance 
metric. The division between the BTB and control animal groups was supported by a bootstrap value of 100%. 



was supported by a 100% bootstrap value after 1,000 
permutations. 

Genome-wide transcriptional profiles generated from 
the PBL of eight M. ^ov/5-infected and eight control ani- 
mals were compared to assess differential gene 



expression between the two sample groups. The micro- 
array analysis revealed a total of 2,960 transcripts, repre- 
senting 2,757 unique genes, as differentially expressed 
between the M. ^ov/5-infected and non-infected control 
animals (adjusted P value < 0.05). In summary, 1,397 
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transcripts (representing 1,281 unique genes) displayed 
increased expression and 1,563 transcripts (representing 
1,476 unique genes) displayed decreased expression in 
the M. Z?ov/5-infected group relative to the control 
group. 

Among the DE genes with immune-related functions 
were genes encoding proinflammatory cytokines and 
other mediators of the host immune response also dis- 
played increased relative expression in the BTB group. 
These included the cluster of differentiation 83 antigen- 
encoding gene {CD83) [+4.92-fold]; the chemokine (C-C 
motif) ligand 2 gene (CCL2) [+2.85-fold]; the chemokine 
(C-X-C motif) ligand 5 gene (CXCLS) [+3.87-fold]; the 
cytotoxic T-lymphocyte-associated protein 4 gene 
(CTLA4) [+3.20-fold]; the chemokine (C-X-C motif) 
receptor 4 gene {CXCR4) [+2.64-fold]; the interleukin 8 
gene {IL8) [+2. 15 -fold] and the interleukin 1 alpha gene 
(ILIA) [ + 1. 89-fold]. The genes encoding the TNF-a 
(TNF) and IFN-y (IFNG) cytokines, two recognised bio- 
markers of M. bovis infection, were not differentially 
expressed between the two groups. 

Reduced relative expression of host immune-related 
genes was also observed in the M, Z?ov/5-infected group 
relative to the control group. Among these were the 
antimicrobial beta-defensin 10 gene (DEFBIO) [-3.38- 
fold]; the triggering receptor expressed on myeloid cells 
1 gene (TREMl) [-1.70-fold]; and the TYRO protein tyr- 
osine kinase binding protein gene [TYROBP) [-1.38- 
fold]. Other immune genes displaying reduced relative 
gene expression in the BTB group included those 
encoding proinflammatory cytokines such as the inter- 
leukin 15 gene {ILIS) [-1.49-fold]; the interleukin 16 
gene {IL16) [-1.44-fold]; and the interleukin 18 gene 
{IL18) [-1.72-fold]. 

Notably, several genes involved in TLR-mediated sig- 
nalling displayed reduced relative expression in the BTB 
group such as the Toll-like receptor 4 gene {TLR4) 
[+2.41-fold]; the Toll-like receptor 2 gene {TLR2) [-1.45- 
fold]; the TLR adaptor protein myeloid differentiation 
primary response gene (88) [MYD88] (-1.31-fold); the 
interleukin- 1 receptor-associated protein kinase 4 gene 
(IRAK4) [-1.11 -fold] and the mitogen-activated protein 
kinase 13 and 14 genes {MAPK13 [-1.93-fold] and 
MAPK14 [-1.15-fold]). The gene encoding the intracellu- 
lar TLR3 protein {TLR3) also displayed reduced relative 
expression in the BTB animals (-1.55-fold). 

Real time quantitative reverse transcription PGR (qRT- 
PCR) analysis and validation of Affymetrix® GeneChip® 
results 

A panel of 23 immune-related genes, including members 
of the interferon signalling pathway, which was recently 
shown to be the most significantly over-represented 
pathway in human patients with active TB [37], were 



selected for real time qRT-PCR analysis. These were 
CASPh CD83, CTLA4, DEFBIO, IFNB, IFNARl, IFNG, 
IFNGRl, IFNGR2, ILIA, IL8, ILIS, JAKl, KIR3DS1, 
MYD88, PTPN2, STATl, STAT2, TLR3, TLR4, TREMl, 
TYK2 and TYROBP, The immune-related function of 
each gene and the results from these analyses are 
detailed in Table 1. 

Thirteen of the 23 genes analysed using real time 
qRT-PCR {CASPl, DEFBIO, IFNARl, ILIS, KIR3DS1, 
MYD88, PTPN2, STATl, STAT2, TLR3, TREMl, TYK2 
and TYROBP) showed significant decreased expression 
{P < 0.05) and five genes {CD83, CTLA4, IFNGRl, ILIA 
and IL8) displayed significant increased expression {P < 
0.05) in the M.bovis-iniected group relative to the con- 
trol animals. No statistically significant differences in 
expression {P > 0.05) were observed between the two 
sample groups for the remaining five genes assayed 
{IFNBl, IFNG, IFNGR2, M/Ciand TLR4), 

Twenty-two of the 23 genes analysed by real time 
qRT-PCR were represented on the microarray; only 
TYK2 was not represented. Gene expression profiles for 
16/22 (73%) of the genes analysed with real time qRT- 
PCR were concordant with the results from the microar- 
ray analysis. Two genes {IFNGR2 and TLR4) displayed 
significant expression differences between the two sam- 
ple groups based on the microarray results, but were 
not significantly different based on real time qRT-PCR 
data analysis. Three genes {IFNARl, IFNGRl and 
PTPN2) displayed significant differences between the 
two groups based on real time qRT-PCR results; how- 
ever, these genes were not differentially expressed 
according to the microarray results. Only one gene 
{STATl) displayed directionally discordant gene expres- 
sion profiles between the two methods. In the M, bovis- 
infected group, STATl showed significant reduced rela- 
tive expression (-1.27-fold) based on real time qRT-PCR 
results but displayed significant increased relative 
expression ( + 1. 28-fold) according to the microarray 
results. The observed discrepancies between the micro- 
array and real time qRT-PCR data may reflect differ- 
ences in the sensitivity of the two analytical methods 
used and/or differences in the mRNA transcripts tar- 
geted by the probes (microarray) and primer pairs (real 
time qRT-PCR) used for the two forms of gene expres- 
sion analysis. 

Analysis of differential gene expression using Ingenuity® 
Systems Pathway Analysis (IPA) 

The total number of DE genes that could be mapped to 
molecules in the Ingenuity® Knowledge Base was 1,869 
from a total of 2,960 DE transcripts. IPA was used to 
categorise these 1,869 DE genes based on their func- 
tional annotation and to assess if a functional gene cate- 
gory contained an over-representation of genes relative 
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Table 1 Gene expression fold-changes between /W.boWs-lnfected (n = 8) and control (n = 8) based on mlcroarray and 
real time qRT-PCR analyses 



Gene 
symbol 


Gene name 


Gene description 


Mean /W.bows-infected 

vs control group 
expression fold-change 
(microarray) 


Mean /W. bows-infected vs 
control group expression 
fold-change (real time 
qRT-PCR) 


Real 
time 
qRT-PCR 
P-value 


CASPl 


Caspase 1 gene 


A member of the cysteine-aspartic acid 
protease. Plays a role in the cell apoptosis 


-1.45 


-1.98 


0.001 


CD83 


Cluster of differentiation 
83 gene 


A cell surface protein found in antigen 
presenting cells. Believed to play a role in 
antigen presentation or the cellular 
interactions that follow lymphocyte 
activation 


+4.92 


+5.88 


< 0.001 


CTLA4 


Cytotoxic T-lympliocyte- 
associated protein 4 
gene 


A member of the immunoglobulin 
superfamily and encodes a protein which 
transmits an inhibitory signal to T cells 


+3.20 


+3.46 


0.001 


DEFBW 


Beta-defensin 10 gene 


Host defence response to bacterial infection; 
has antimicrobial activity 


-3.38 


-12.21 


0.018 


IFNBl 


interferon, beta 1, 
fibroblast 


Cytokine activity 


Not DE 


Not DE 


> 0.05 


IFNARl 


Interferon (alpha, beta 
and omega) receptor 1 
gene 


A type 1 membrane protein that forms one 
of the two chains of a receptor for 
interferons alpha and beta 


Not DE 


-2.20 


< 0.001 


IFNG 


Interferon gamma gene 


A soluble cytokine with antiviral, 
immunoregulatory and anti-tumor properties. 
A potent activator of macrophages 


Not DE 


Not DE 


> 0.05 


IFNGRl 


Interferon gamma 
receptor 1 gene 


Forms a heterodimer with interferon gamma 
receptor 2. Involved in binding of interferon 
gamma 


Not DE 


+2.23 


< 0.001 


IFNGR2 


Interferon gamma 
receptor 2 gene 


Forms a heterodimer with interferon gamma 
receptor 1. Involved in binding of interferon 
gamma 


-1.32 


Not DE 


> 0.05 


ILIA 


Interleukin 1, alpha 
gene 


A member of the interleukin 1 cytokine 
family produced by macrophages. Involved 
in various immune responses, inflammatory 
processes, and haematopoiesis 


+1.89 


+2.47 


0.009 


IL8 


Interleukin 8 gene 


A chemokine that mediates the 
inflammatory response 


+2.15 


+4.01 


< 0.002 


IL15 


Interleukin 15 gene 


A cytokine that regulates T and natural killer 
cell activation and proliferation 


-1.49 


-2.03 


< 0.001 


JAKl 


Janus kinase 1 gene 


A widely expressed membrane-associated 
phosphoprotein involved in the interferon- 
alpha/beta and -gamma signal transduction 
pathways 


Not DE 


Not DE 


> 0.05 


KIR3DS1 


Killer cell 

immunoglobulin-like 
receptor, three domains, 
short cytoplasmic tail, 1 


A transmembrane glycoprotein expressed by 
natural killer cells and some T cells. Involved 
in regulation of the immune response 


-1.44 


-2.63 


0.005 


MYD88 


myeloid differentiation 
primary response gene 
(88) 


A cytosolic adapter protein that functions as 
a signal transducer in the interleukin-1 and 

~rlll'l j_ • 11* j_l 

Toll-like receptor signalling pathways 


-1.31 


-2.55 


< 0.001 


PTPN2 


protein tyrosine 
phosphatase, non- 
receptor type 2 gene 


A member of the protein tyrosine 
phosphatase (PTP) family that is involved in 
a variety of cellular processes including cell 
growth, and differentiation 


Not DE 


-1.43 


< 0.001 


STATl 


signal transducer and 
activator of transcription 
1, 91kDa gene 


A transcriptional activator protein activated 
in response to cytokines and growth factors 


+1.28 


-1.27 


< 0.001 


STAT2 


signal transducer and 
activator of transcription 
2, llSkDa gene 


A transcriptional activator protein activated 
in response to cytokines and growth factors 


-1.23 


-1.86 


< 0.001 
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Table 1 Gene expression fold-changes between /W.bows-infected (n = 8) and control (n = 8) based on mlcroarray and 
real time qRT-PCR analyses (Continued) 



TLR3 Toll-like receptor 3 gene An intracellular pathogen recognition -1.55 -1.63 0.001 

receptor (PRR) that largely recognises viral 
pathogen-associated molecular patterns 
(PAMPs). Mediates the production of 
cytokines necessary for the development of 
effective immunity 

TLR4 Toll-like receptor 4 gene A cell-surface PRR that recognises PAMPs +2.41 Not DE > 0.05 

expressed by infectious agents. Mediates the 
production of cytokines necessary for the 
development of effective immunity 

A receptor expressed on myeloid cells upon -1.70 -2.81 0.002 

microbial infection. Amplifies neutrophil and 
monocyte-mediated inflammatory responses 
stimulating release of pro-inflammatory 
chemokines and cytokines 



TYK2 Tyrosine kinase 2 gene 


A member of the tyrosine kinase. A 
component of both interferon signalling 
pathways 


Not represented on array -l.( 


50 < 0.001 


TYROBP JYRO protein tyrosine 


A transmembrane protein involved in cell 


-1.38 -1.^ 


58 < 0.001 


kinase binding protein 


signalling 






gene 









TREMl Triggering receptor 

expressed on myeloid 
cells 1 gene 



Gene descriptions taken from GeneCards version 3 [82] and the NCBI Entrez Gene database [83]. For the microarray data, genes with fold-change differences In 
expression were significant after adjustment for multiple testing using the Benjamlnl-Hochberg method [31] (adjusted P < 0.05). 'Not DE' Indicates a gene was 
not differentially expressed between the two sample groups. 



to the microarray reference gene list. This analysis 
showed that the top functional category observed for the 
DE genes used for IPA was the inflammatory response, 
which contained 241 genes with P-values ranging from 
1.66 X 10"^^ to 1.64 X 10'^ (Table 2). The inflammatory 
response category was further sub-divided with the affects 
immune response subcategory containing the most mole- 
cules (138 DE genes) [Figure 3]. Further inspection of the 
individual genes within this subcategory revealed that 
there was an over-representation of genes displaying 
reduced relative expression (89 genes) compared to DE 
genes showing increased relative expression in the BTB 
animals (49 genes). This observation was in contrast to 
the other top functional categories where similar num- 
bers of genes showing increased and decreased relative 
expression were reported. 

Canonical molecular pathways associated with M, 
^ov/5-infection were analysed using IPA. These canoni- 
cal pathways were ranked according to P-value, which 
represents the significance of the association between a 
specific pathway and the genes in the input data set. 
The majority of the top ranking IPA-identified canonical 
pathways were involved in cell signalling and communi- 
cation associated with host innate and adaptive immune 
responses. Natural killer cell signalling and communica- 
tion between innate and adaptive immune cells were 
identified as the top ranking canonical pathways. In 
addition, TREMl signalling dendritic cell maturation, 
JAK-ST AT signalling, T cell signalling, IL6 signalling, 
chemokine signalling and TLR signalling were among the 



top twenty IPA-identified canonical pathways (Table 3). 
Based on the well-documented role of TLR signalling in 
mycobacterial infection [13,38-42], this canonical path- 
way overlaid with gene expression results is shown in 
Figure 4. 

Discussion 

The implementation of surveillance and management 
programmes has done much to reduce the incidence 
and prevalence of BTB over the past number of decades; 
however, M. bovis infection remains an important live- 
stock disease worldwide. This is due, in part, to well- 
documented limitations of the currently available diag- 
nostics tests (such as the SICTT and IFN-y tests) lead- 
ing to a failure to detect all infected animals [43,44]. 

In recent years, research has shifted from a focus on 
protein-based diagnostics to functional genomics tech- 
nologies that interrogate the host transcriptome in 
response to M, bovis infection. In particular, microarray 
technologies coupled with the rapid development of 
more sophisticated bovine genome resources has 
enabled high-resolution analyses of the genes and cellu- 
lar pathways governing the host response to infection 
with M. bovis [23,25,45-47]. In the present study, we 
have compared the transcriptomes of PBL from non- 
infected control animals with actively-infected BTB ani- 
mals using a high-density genome-wide bovine microar- 
ray platform. 

Modulation of the host PBL transcriptome in response 
to M. bovis infection was evident from the large number 
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Table 2 Gene ontology (GO) categories Identified using IPA 



Gene ontology category 


P-value range 


Number of genes 


Inflammatory response 


1 .66 X 1 0 - 1 .64 X 1 0 


241 


Cellular development 


8.95 X 1 0 - 1 .57 X 1 0 


308 


Cellular growth and proliferation 


1.65 X 10 - 1.57 X 10 


333 


Haematological system development and function 


1 .65 X 1 0 - 1 .64 X 1 0 


279 


Haematopoiesis 


2.60 X 1 0 - 1 .57 X 1 0 


1 73 


Tissue morphology 


1 .78 X 1 0 -1.19x10 


1 14 


Cellular function and maintenance 


5.44 X 1 0 - 1 .65 Xl 0 


1 38 


Cell death 


9.57 X 10 - 1.60 X 10 


419 


Cell-mediated immune response 


1 .06 X 1 0 - 1 .25 X 1 0 


1 09 


Connective tissue disorders 


1 .72 X 1 0 - 5.94 X 1 0 


234 


Immunological disease 


1 .72 X 1 0 - 1 .57 X 1 0 


344 


Inflammatory disease 


1 -7") 1 r\"6 1 -10 ^y 1 r\"2 

1 .72 X 1 0 - 1 .28 X 1 0 


364 


Skeletal and muscular disorders 


1 .72 X 1 0 - 1 .22 X 1 0 


235 


Infectious disease 


1 .82 X 1 0 -1.17x10 


223 


Cell-to-cell signalling and interaction 


8.33 X 1 0 - 1 .64 X 1 0 


210 


Dermatological diseases and conditions 


1 .04 X 1 0 - 1 .57 X 1 0 


1 07 


Free radical scavenging 


1.23 X 10 - 1.23 X 10 


14 


Molecular transport 


1.23 X 10 - 7.43 X 10 


86 


Cell signalling 


2.40 X 1 0 - 7.21 X 1 0 


1 33 


Respiratory disease 


2.86 X 1 0 - 1 .24 X 1 0 


1 09 


Cellular compromise 


r\r\ \ y 1 r\~5 1 r ~7 ^ , 1 /^~2 

3.09 X 1 0 - 1 .57 X 1 0 


41 


Genetic disorder 


3.39 X 1 0 - 1 .64 X 1 0 


1 04 


Immune cell trafficking 


3.51 X 1 0 - 1 .64 X 1 0 


163 


Humoral immune response 


5.08 X 1 0 -1.19x10 


78 


Organismal injury and abnormalities 


9.51 X 10 - 1.05 X 10 


81 


Tissue development 


1.19x10 - 1 .05 X 1 0 


1 08 


Haematological disease 


1 -70 ^y 1 r\"4 1 T r 1 r\"2 

1 .78 X 1 0 - 1 .35 X 1 0 


162 


Gastrointestinal disease 


1 .88 X 1 0 - 5.51 X 1 0 


42 


Hepatic system disease 


1 .88 X 1 0 - 1 .64 X 1 0 


89 


Lymphoid tissue structure and development 


1 .94 X 1 0 - 5.94 X 1 0 


5 


Antigen presentation 


~> A ^ ^ y 1 r\~4 r A ^ y 1 r\~2 

3.42 X 1 0 -1 .64 X 1 0 


83 


Cellular movement 


3.58 X 1 0 - 1 .63 X 1 0 


1 57 


Vitamin and mineral metabolism 


3.59 X 1 0 - 7.21 X 1 0 


74 


Antimicrobial response 


3.92 X 10 - 1.28 X 10 


32 


Hypersensitivity response 


5.86 X 1 0 - 8.60 X 1 0 


31 


Infection mechanism 


r\ c>r\ ^y 1 r\"4 -[ c~i ^y i r\~2 

9.89 X 1 0 - 1 .57 X 1 0 


39 


Gene expression 


1 .07 X 1 0 - 1 .57 X 1 0 


254 


Hair and skin development and function 


1.23 X 10 - 1.57 X 10 


22 


Post-translational modification 


n n-7 K ^ 1 r\"3 n -70 K ^ 1 r\"3 

2.27 X 10 - 9.78 X 10 


84 


Cancer 


TCI \y 1 r>"3 1 c-7 \y 1 r\"2 

2.5 1 X 1 0 - L5/ X 1 0 


385 


Endocrine system development and function 


2.65 X 10"^ - 2.65 x 10"^ 


21 


Carbohydrate metabolism 


3.77 X 10"^ - 3.77 x 10"^ 


6 


Cell morphology 


3.77 X 10"^ - 1.28 X 10"2 


23 


Lipid metabolism 


3.77 X 10"^ - 7.43 x 10"^ 


17 


Neurological disease 


3.77 X 10"^ - 1.47 X 10"^ 


49 


Small molecule biochemistry 


3.77 X 10"^ - 7.43 x 10"^ 


23 


Metabolic disease 


3.91 X 10"^ - 3.91 X 10"^ 


7 
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Table 2 Gene ontology (GO) categories Identified using IPA (Continued) 



Embryonic development 


4.59 X 10"^ - 1.57 X 10"^ 


21 


Nervous system development and function 


4.59 X 10"^ - 1.18 X 10"^ 


8 


Reproductive system disease 


5.84 X 10"^ - 7.54 x 10"^ 


79 


Cell cycle 


5.94 X 10"^ - 1.18 X 10"^ 


20 


Endocrine system disorders 


5.94 X 10"^ - 5.94 x 10"^ 


3 


Opthalmic disease 


5.94 X 10"^ - 6.60 X 10"^ 


7 


RNA post-transcriptional modification 


6.04 X 10"^ - 6.04 X 10"^ 


13 


Cardiovascular system development and function 


6.68 X 10"^ - 1.30 xlO"^ 


19 


Cell assembly and modification 


6.68 X 10"^ - 1.28 X 10"^ 


14 


Organ development 


6.68 X 10"^ - 1.19 X 10"^ 


6 


Renal and urological system development and function 


1.18 X 10"^ - 1.30 X 10"^ 


11 


Reproductive system development and function 


1.18 X 10"^ - 1.18 X 10"^ 


4 



The top ranking GO categories identified by IPA are listed according to P-values. P-value ranges are based on subcategory P-values within each parent term. 

of DE genes between the two experimental groups. Sta- the M. Z?ov/5-infected and control animals (Figure 1) 

tistical analyses of the microarray data identified a total may have contributed to the gene expression changes 

of 2,757 DE genes. Of these, 1,281 (46%) showed detected between the two experimental groups. Also, 

increased expression and 1,476 (54%) displayed the haematology analyser results only provided a general 

decreased expression in the BTB group compared to the description of the PBL cell subsets and do not provide 

control animals. It is important to note, however, that information concerning T lymphocyte subsets in the 

the differences in cell subpopulations observed between infected and control animal groups. In addition, the cell 



□ Increased □ Decreased 




Figure 3 Sub-categories of the top ranlcing IRA-identified inflammatory response gene ontology (GO) category. The numbers of genes 
displaying increased and decreased relative expression for affects immune response, the top ranking subcategory, are shown. The number of 
differentially expressed genes within each functional subcategory is indicated. 
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Table 3 Top-ranking canonical pathways identified using IPA 



Canonical Pathway Name 


P-value 


Ratio 


Natural killer cell signalling 


3.06 X 1 0'^ 


25/ \ 1 Z (0.20b) 


Communication between innate and adaptive immune cells 


4.1 1x10'^ 


1 7/89 (0.1 91 ) 


TREMl signalling 


2.20 X 1 0'^ 


15/69 (0.217) 


Dendritic cell maturation 


4.96 X 1 0 


29/1 74 (0.167) 


Cysteine metabolism 


7.43 X 10'^ 


9/90 (0.100) 


1 A 1/ /c — TAT ^; ^ 1 1 ; ^ 

JAK/blAI Signalling 


8.22 X 10'^ 


15/64 (0.234) 


NRF2-mediated oxidative stress response 


0.98 X 1 0 


35/1 83 (0.1 91 ) 


T cell receptor signalling 


1 .08 X 10 


20/107 (0.187) 


IL6 signalling 


1 .52 X 1 0'^ 


20/93 (0.215) 


CCR5 signalling in macrophages 


1 .84 X 10 


16/92 (0.1 74) 


Chemokine signalling 


1 .92 X 1 0'^ 


1 r /-7r ir\ ^r\r\\ 

15/75 (0.200) 


Calcium-induced T-lymphocyte apoptosis 


2.07 X 10'^ 


12/66 (0.182) 


IL-1 7 signalling 


2.21 X 10'^ 


16/74 (0.216) 


Prolactin signalling 


2.21 X 10'^ 


16/75 (0.213) 


Synaptic long term potentiation 


2.22 X 10'^ 


1 n /I IT ir\ 1 r o\ 

19/113 (0.1 68) 


Toll-like receptor signalling 


2.63 X 10'^ 


1 1/54 (0.204) 


FLT3 signalling in hematopoietic progenitor cells 


2.77 X 10'^ 


15/74 (0.203) 


Systemic lupus erythematosus signalling 


2.81 X 10'^ 


20/151 (0.132) 


Renin angiotensin signalling 


3.01 X 10'^ 


21/1 20 (0.1 75) 


Phospholipase C signalling 


3.28 X 1 0'^ 


39/253 (0.154) 


Oncostatin M signalling 


3.36 X 1 0'^ 


9/35 (0.257) 


Thyroid cancer signalling 


3.36 X 1 0'^ 


9/42 (0.214) 


B cell receptor signalling 


3.72 X 10'^ 


25/154 (0.162) 


1 r\ "horror/^ r\ ci/^r\::3 

IIILcllclUII biyiidiiiiiy 


Q Q7 V 1 n"^ 

o.y / A 1 U 




Production of nitric oxide and reactive oxygen species in macrophages 


4.13 X 10"^ 


28/185 (0.151) 


NF-kB signalling 


4.17 X 10"^ 


25/152 (0.164) 


Notch signalling 


4.29 X 10"^ 


8/43 (0.186) 


IL10 signalling 


4.36 X 10"^ 


13/70 (0.186) 


P38 MARK signalling 


4.57 X 10"^ 


17/97 (0.175) 


Role of NFAT in regulation of the immune response 


4.70 X 10"^ 


29/196 (0.148) 


Parkinson's signalling 


5.00 X 10"^ 


5/1 7 (0.294) 



The top ranking canonical pathways identified by IPA are listed according to P-values. The ratio indicates the number of differentially expressed genes involved 
in each canonical pathway divided by the total number of genes within each pathway as per the IPA Knowledge Base. 



subset results presented here differ from previous work 
performed by us [25], most likely due to the different 
cell sample types examined. 

Analysis of the DE genes using IPA provided informa- 
tion regarding the immunobiology of active BTB. The 
highest ranking functional category identified using IPA 
was inflammatory response and the highest ranking sub- 
category within this category-<^^ec^5 immune response- 
revealed a marked bias in the number of genes display- 
ing a decrease in relative expression (64.5%) compared 
to those showing an increase in relative expression 
(35.5%) in PBL from the BTB group. 

Previous work by our research group demonstrated 
that suppression of host gene expression was associated 



with active M. bovis infection in cattle [25]. This earlier 
work involved the comparison of RNA isolated from 
PBMC of M. ^ov/5-infected and control animals using 
an immuno-specific bovine cDNA microarray (BOTL-5). 
The results presented here, based on the analysis of 
PBL-derived RNA using a genome-wide microarray, 
lend further support to our previous study. Indeed, pub- 
lished investigations by other workers suggest that tran- 
scriptional suppression is a common feature of 
mycobacterial infection in mammals [48-50]. 

Further inspection of the inflammatory response func- 
tional category in IPA identified several genes that were 
previously reported to be differentially expressed in cat- 
tle and other mammalian species infected with 
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Figure 4 Differential gene expression in the TLR signalling pathway. Genes within tine TLR signalling pathway showing differential 
expression are highlighted in colour. Colour intensity indicates the degree of increased (red) or decreased (green) relative expression in the BTB 
group compared to the control animal group. A white colour indicates genes that were not differentially expressed and entities coloured grey 
represent microbial PAMPs. 



mycobacterial pathogens [23,40,42,51,52]. For example, 
microarray analysis showed that TLR2 and TLR4 dis- 
played contrasting expression patterns between PBL 
from the two groups: TLR2 showed decreased expres- 
sion and TLR4 showed increased expression relative to 
the control animal group. We have previously observed 
decreased expression of TLR2 in PBMC from actively 
infected BTB cattle using the immuno-specific BOTL-5 
cDNA microarray; however, contrary to the results of 
the present study, TLR4 also showed decreased 



expression with the BOTL-5 cDNA microarray in 
actively infected animals [25]. 

The gene expression results obtained by Meade and 
colleagues using PBMC from M. ^ov/5-infected and con- 
trol non-infected animals were also used to identify a 
panel of 15 genes predictive of disease status [25]. Four 
of these genes were found to be similarly differentially 
expressed in the current study: UNC84B (now SUN2), 
GAK SFPQ and NRPl. Four other of the 15 genes iden- 
tified previously {TBKl 28S [now RN28S1], GPR98 and 
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an anonymous BOTL clone [BOTL0100013_F01]) were 
not present on the Affymetrix® GeneChip® Bovine Gen- 
ome Array. However, the seven remaining genes 
[NCORl, PPP2RSB, UCP2, ZDHHC19, NFKBl, NRM 
and FGFRl) were not differentially expressed in the PEL 
samples from M. Z?ov/5-infected and control non-infected 
animals used for the present study. This discordance 
may be due to a number of factors, including: the blood 
cell sample types used (PEL versus PEMC); differences 
in sensitivity between the two types of microarray (the 
single-colour in 5/^r/-synthesised Affymetrix® GeneChip® 
versus a dual-colour spotted cDNA array [53,54]); and 
the requirement for more stringent control of the FDR 
with larger numbers of genes (24,072 probe sets versus 
1,391 duplicate spot features). 

The role of TLR molecules in the recognition of 
mycobacterial PAMPs is well established 
[11,38,39,41,42,55-57]. TLR2 and TLR4 activation sig- 
nals are Unked to the interleukin-1 receptor-associated 
protein kinases (IRAKs) through the adaptor molecule, 
myeloid differentiation primary response protein 88 
(MYD88), which triggers a downstream protein signal- 
ling cascade involving tumour necrosis factor receptor- 
associated factor 6 (TRAF6) and mitogen-activated pro- 
tein kinases (MAPKs) [58,59]. This cascade culminates 
in the expression of many NF-^E-inducible genes, 
including CCL2, ILIB, IL12, IL18 and TNF, causing nat- 
ural killer (NK) and T cells to release IFN-y and TNF-a, 
which ultimately results in granuloma formation [60]. 

In the present study, several TLR-mediated proinflam- 
matory cytokines and signalling molecules were differen- 
tially expressed in the ETE group compared to the non- 
infected control animals. These included CCL2 
(increased), CXCR4 (increased), CXCL5 (increased), 
ILIA (increased), IL8 (increased), IL18 (decreased), 
IRAK4 (decreased), MAPK6 (increased), MAPK13 
(decreased), MAPK14 (decreased) and MYD88 
(decreased). This was also supported by canonical path- 
way analysis using IPA, which identified TLR signalling 
as a molecular pathway affected by M. bovis infection. 

These results suggest that genes encoding TLR- 
mediated signalling pathway molecules have a role in 
governing the host response to ETE and may also serve 
as targets for immuno-subversion by M. bovis. For 
example, genes encoding several innate immune recep- 
tors and chemokines (such as TLR4, CD83, CCL2, 
CXCR4, CXCLS, ILIA and 718)— several of which parti- 
cipate in the initiation of a T cell response during infec- 
tion [61-64]— showed increased relative expression in 
the ETE animal group. Transcriptional profiles suggest- 
ing initiation of a T cell response are supported by the 
comparative analysis of the PEL cell populations in the 
two animal groups; a significant increase in the mean 
number of lymphocytes and a significant decrease in the 



mean number of monocytes were observed in the ETE 
group relative to the control animals. This difference in 
the PEL cell composition may represent recruitment of 
host cytotoxic lymphocytes for the destruction of 
infected monocytes in the control of M. bovis infection 
[3,65]. 

It is important to note, however, that the observed 
decreased expression of host PRR genes (such as TLR2) 
and the genes encoding their associated adaptor and sig- 
nalling pathway molecules (such as MYD88, IRAK4, 
MAPK13 and MAPK14) may indicate that the adaptive 
response in ETE animals is inferior due to the repres- 
sion of these innate immune genes. Indeed, previous 
work has proposed that mycobacterial antigens, such as 
the early secreted antigenic target protein 6 (ESAT-6) 
protein, attenuates the host innate immune response by 
inhibiting MYD88-IRAK4 binding, thus causing suppres- 
sion of NF-A^E-induced transcription of upstream genes 
required for T cell response initiation [66]. These work- 
ers also demonstrated that activation of v-akt murine 
thymoma viral oncogene homolog kinases (AKTs) is 
necessary to prevent MYD88-IRAK4 complex formation. 
Notably, the AKT2 gene displayed increased relative 
expression (+ 1.22-fold) in the ETE animal group in the 
present study. 

Repression of host innate immune genes that elicit an 
adaptive response to M. bovis infection is further sup- 
ported by the analysis of genes belonging to the inter- 
feron signalling pathway, which has been shown to have 
a role in human tuberculosis [37,67-71]. IFN-y is 
secreted by NK cells and CD4"^ T cells upon activation 
by IL-12 produced by infected macrophages. IFN-y 
recruits additional macrophages to the site of infection 
while also providing the stimulus for activating microbi- 
cidal functions in infected macrophages [14,71-73]. IFN- 
y also induces MHC class II gene expression in infected 
macrophages by signalling through its receptor (IFN-y- 
receptor) [74-76]. This stimulates the JAK-STAT path- 
way, resulting in induction of transcriptional activators 
of MHC class genes, such as the MHC class II transacti- 
vator gene {CUT A), Mycobacterial antigen presentation 
via MHC class II molecules is critical for the recruit- 
ment of additional CD4^ T cells and the formation and 
maintenance of granulomas [38]. 

The results from the current study support a role for 
interferon signalling pathways during M. bovis infection. 
The genes encoding interferon (alpha, beta and omega) 
receptor 2 {IFNAR2); interferon gamma receptor 2 
{IFNGR2); interferon-induced protein with tetratricopep- 
tide repeats 2 {IFIT2); interferon-induced protein with 
tetratricopeptide repeats 5 {IFIT5); interferon-induced 
transmembrane protein 3 {IFITM3); protein tyrosine 
phosphatase, non-receptor type 2 {PTPN2); and signal 
transducer and activator of transcription 2 {STAT2) 
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displayed differential expression in the BTB animals 
based on the microarray and/or real time qRT-PCR 
analyses. 

These findings suggest that, in addition to the target- 
ing of TLR-mediated signalling pathways, M. bovis may 
also target genes involved in the IFN-signalling pathway, 
resulting in an attenuated T cell response that enables 
mycobacterial survival and disease progression. It is 
tempting to speculate that suppression of IFN-signalling 
in response to M. bovis infection may result in the 
impairment of the antigen presenting process required 
for adaptive immunity; however, further work is 
required to investigate this possibility. Notably, the gene 
encoding IFN-y (IFNG) was not differentially expressed 
in the current study, despite the BTB animals having 
tested positive for increased IFN-y based on the Bovi- 
GAM® assay. However, it is important to note that 
unlike the blood samples used for the BoviGAM® assay 
in the current work, the PBL fraction from which the 
RNA was derived in this study was not stimulated with 
protein purified derivative of tuberculin (PPD), which is 
required to elicit IFN-y secretion [4]. In addition, con- 
trary to previous results obtained by Meade and collea- 
gues [25] we did not detect differential expression of the 
TNF gene between the M. ^ov/5-infected and control 
animals examined here. The most likely explanation for 
this apparent discrepancy is the different cell sample 
types used for gene expression analyses (PBL versus 
PBMC). 

IPA canonical pathway analysis identified a number of 
DE genes which, to our knowledge, have not previously 
been reported to be involved in the host response to 
tuberculosis in cattle or other mammalian species. 
These included CTLA4 and TLR3, TLR3 encodes an 
intracellular PRR involved in the recognition of viral- 
derived nucleic acids [77]. In the present study, reduced 
relative TLR3 expression in BTB animals (-1.55-fold) 
may suggest some hitherto unknown role for this PRR 
in intracellular mycobacterial infection. In support, we 
have observed significant differential expression of TLR3 
in bovine monocyte-derived macrophages (MDM) sti- 
mulated in vitro with M. bovis when compared to non- 
stimulated control MDM (unpubUshed data). CTLA4 
encodes an inhibitor of the T cell-mediated response 
[78,79] and this gene displayed increased relative expres- 
sion (-1-3.20-fold, P < 0.05) in the M. bovis infected ani- 
mals in the present study. The observed increased 
relative expression of CTLA4 may reflect a mechanism 
of immuno-modulation used by M. bovis to subvert a 
host T cell response. 

Finally, hierarchical clustering analysis was performed 
here using a total of 5,388 genes that passed the infor- 
mative probe filtering criteria. This analysis unambigu- 
ously differentiated animals on the basis of their disease 



status. This result suggests that genome-wide expression 
profiling of PBL from BTB animals can be used to 
enable the identification of suitable transcriptional mar- 
kers for the detection of infected animals within herds 
and augment current surveillance strategies in countries 
where control programmes have been implemented 
[21,22]. However, further work using PBL samples from 
additional animals infected with M. bovis and other 
microbial pathogens will be required to identify and 
validate robust M. ^ov/5-specific transcriptional signa- 
tures of infection. 

Conclusions 

The results presented here support the hypothesis that 
repression of immune-related genes is an important fea- 
ture of mycobacterial infections [25,48-50,80,81]. In par- 
ticular, the gene expression results obtained suggest that 
M. bovis infection may target the innate immune cellu- 
lar pathways necessary for the initiation of the appropri- 
ate T cell response. Notably, analysis of the cell 
populations present in the PBL from the BTB animals 
showed an increase in the number of lymphocytes rela- 
tive to the control animals, suggesting that the actively- 
infected BTB animals do mount a T cell response. How- 
ever, it is possible that the T cell response elicited by 
these animals is compromised, resulting in disease pro- 
gression. Indeed, failure of the adaptive immune 
response to contain the mycobacterial infection is 
regarded as the primary cause of the development of 
active tuberculosis from a latent state of infection [74]. 
Finally, cluster analysis using all informative mRNA 
transcripts permitted a clear delineation between healthy 
and infected animals. These results demonstrate that 
functional genomics approaches based on transcriptional 
profiling can be used to supplement current protein- 
based diagnostics for BTB. 
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Additional file 1: Perl script used for bootstrapping of cluster 
analysis results in Figure 2. 

Additional file 2: Table S1: Oligonucleotide primers used for real 
time qRT-PCR validation of microarray results. 
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